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I. INTRODUCTION 

The nature of the magnetic excitation spectrum of the 
cuprates is an important topic in its own right. But it 
has many imphcations beyond this as well. It is of spe- 
cial relevance to those theories that propose a 'magnetic' 
origin for cuprate superconductivitj^. And it is of par- 
ticular interest to the current debate whether magnetic 
excitations or phonons are responsible for certain strong 
coupling features observed in tunnelling, infrared conduc- 
tivity, and angle resolved photoemission spectr 

Although there are a large number of theories for ex- 
plaining the observed magnetic excitation spectrum, to a 
first approximation, these theories can be collapsed into 
two groups. The first are linear response calculations 
based on a two dimensional Fermi surface and the pres- 
ence of a d-wave energy gap. In this group are included 
calculations which go beyond RPA (such as FLEX^) and 
the generalization to the particle-particle channel of the 
SO (5) approach^. The other group are those based on 
coupled spin ladders^-^. They assume phase segrega- 
tion of the material into undoped antiferromagnetic do- 
mains separated by one-dimensional 'stripes' containing 
the doped holes. They have the advantage of treating the 
full quantum mechanical nature of the spins, but have 
the disadvantage that fermionic degrees of freedom are 
neglected. 

It is not the intent here to discuss the relative merits 
of one approach versus the other. In some sense, both 
theories are different limits of a more complex theory 
which properly treats spin and charge degrees of freedom. 
Rather, I wish to discuss what success linear response 
calculations have in regards to describing magnetic exci- 
tation properties of the cuprates, but as well discuss their 
limitations. In particular, I wish to address the question 
of whether there is a 'universal' behavior of the spin exci- 
tation spectrum based on the RPA results. Ultimately, I 
hope these results will have some relevance to a final un- 
derstanding of the magnetic behavior of these fascinating 
materials. 



II. METHODOLOGY 



The RPA expression for the interacting susceptibility 



IS 



Xo(g,^) 

1 - Uxo{q,i^) 
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In this equation, xq is the polarization bubble con- 
structed from bare Greens functions, and U is the effec- 
tive Hubbard interaction which results from projection 
onto the single band subspace (that is, the band that is 
the antibonding combination of copper (1^2 _y2 and oxy- 
gen Pj, and Py orbitalg^). 

It is often stated that RPA is a weak coupling ap- 
proach, applicable only to heavily overdoped materials 
which exhibit Fermi liquid like behavior-'^. But it was 
shown by Schrieffer, Wen, and Zhang^^ that RPA repro- 
duces the superexchange J and the resulting spin wave 
dispersion of the undoped material. Obviously, there are 
important quantum corrections to these results, but these 
do not impact the overall correctness of the approach. 
One imp ortant point to remark is that the calculation of 
Ref. [ll| was self-consistent in the sense that the gapless 
nature of the spin excitations at q = (tt, tt) is a direct 
consequence of the mean-field equation for the Hubbard 
gap. Such self-consistency is difficult to implement in the 
doped case because of the added complication of having 
to calculate the screening caused by the doped holes. 

There are two ways to proceed for the doped case. 
Continue to use bare Greens functions and replace U by 
an effective Ueff which accounts for screenin g^^i^'^ . Or 
take the infinite U limit, replace the Greens functions 
by those of the t- J model, and replace U by J^f / where 
Jeff takes into account the reduction of J by the intro- 
duction of doped holes (in the t-J approach, one is forced 
to scale J down to prevent a magnetic instability from 
occurring^^). An important additional remark is that the 
screened interaction is now dependent on transfered mo- 
mentum, q. This is particularly clear in the t-J model, 
where J{q) = —J{cos{qxa) + cos{qya))/2 due to exchange 
between near neighbor copper sites. 
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One limitation of the RPA approach is the use of bare 
Greens functions as opposed to dressed ones. But it is 
well known that using dressed Greens functions means 
that vertex corrections must be includedi^. The neglect 
of such vertex corrections typically leads to worse results 
than the RPA. 

The approach taken in this paper will be somewhat 
more phenomenological. The bare Greens functions will 
be based on dispersions derived from angle resolved pho- 
toemission. The resulting Ueff will be taken as an ad- 
justable constant, designed to reproduce the location of 
the S=l resonance mode at q = (tt, tt). Calculations have 
been performed as well assuming Uef / has the same q de- 
pendence as J{q). In a more complete theory, one might 
work instead within the context of dynamical mean field 
theory and its cluster generalization^^. 

The heart of the RPA calculation turns out to be the 
structure of the bare polarization bubble, xo- For a su- 
perconductor this isi^ 

El 1 Q efc£fc+g + AfcAfc+g f{Ek+q)-f{Ek) 
^,2^ EkEk+q ^ uj - (Ek+q - Ek) + iS 

1 _ ek^k+q + AkAk+q 1 - f{Ek + q) - fjEk) 

4 ^ EkEk+q ' LJ + {Ek+q + Ek) + i6 

1 _ euek+q + AfcAfc+, . f{Ek+q) + f{Ek) ~ 1 

4 ^ EkEk + q ' LO - [Ek+q +Ek)+i5^ ' 

where Ek = \/ ^1 + Af.. The coherence factors play a 
critical role in the results^i. For an s-wave supercon- 
ductor, the coherence factor (in term 3 of Eq. 2) van- 
ishes on the Fermi surface. The net result is that there 
are no spin collective modes in this case. On the other 
hand, for an order parameter satisfying the condition 
A(fc + q) — —A{k), the coherence factor becomes maxi- 
mal on the Fermi surface (equal to 2). As a consequence, 
the imaginary part of the bubble has a step discontinuity 
at threshold, and thus the real part has a logarithmic di- 
vergence. This divergence guarantees the appearance of 
a collective mode below threshold. This collective mode 
energy at g = (tt, tt) will be denoted as the resonance 
energy. 

There is a common misconception that the appearance 
of the collective mode requires the presence of a d-wave 
superconducting gap. In fact, the only requirement is 
that the coherence factors do not vanish on the Fermi sur- 
face. As an example, imagine that the pseuodgap phase 
is a phase disordered version of a superconductor. To a 
first approximation, we can then set < AkAk+q >=0^. 
In this case, the coherence factor is now unity on the 
Fermi surface independent of the value of fc, and one 
finds a collective mode below threshold as in the case 
of a d-wave superconductor. As the coherence factor is 
now half that of the superconductor, then larger values of 
Ueff are needed to obtain similar collective mode ener- 
gies. In this paper, a d-wave order parameter of the form 
A(fc) = Ao(cos(A:a;a) — cos{kya))/2 is assumed, although 



for some calculations the phase disordered approximation 
was invoked as well. 

Finally, some technical remarks. Eq. 2 is solved by re- 
placing S by some finite F and doing a simple sum over 
the Brillouin zone. The smaller F is, the more k points 
are needed in the sum. Results have been generated for F 
ranging between 0.1 and 2 meV, with zone meshes rang- 
ing from 400 by 400 to 4000 by 4000, but in this pa- 
per, results are shown for r=2 meV (400 by 400 mesh). 
All calculations were performed for a temperature of 1 
meV. All susceptibilities quoted here are states per eV 
per Cu02 formula unit and should be multiplied by the 
matrix element J2a d'^ f^'si'^l^^l'^)'^ ~ ^Ms to compare to 
experiment. 



III. RESULTS - NORMAL STATE 

The most important input into the RPA calculations 
is the assumed form of the fermionic dispersion. For in- 
stance, a well studied case is a tight binding model with 
just near neighbor hopping^. When considering the re- 
sponse at q — (7r,7r), the threshold for the polarization 
bubble corresponds to exciting from the node of the d- 
wave order parameter, kjf, to an unoccupied state, k^+q, 
along the zone diagonal. As such, the threshold is inde- 
pendent of the d-wave energy gap. The Fermi surfaces 
observed by angle resolved photoemission, though, dif- 
fer considerably from those predicted by a near neigh- 
bor tight binding model. In particular, the observed 
Fermi surfaces are characterized by 'hot spots', that is, 
points on the Fermi surface which satisfy the condition 
ffc = Cfe+g = 0. In this case, the threshold of the bubble 
at q = (tt, tt) is determined by twice the energy gap at 
the hot spots. 

When considering the response for a general q, there 
is even more sensitivity to the assumed fermionic disper- 
sion. To illustrate this point, in Fig. 1, I show the real 
part of the normal state (A = 0) bubble at w = for 
four different dispersions along the three symmetry axes 
of the 2D Brillouin zone. The first dispersion (tbl) is that 
based on an early tight binding fit to ARPES spectra on 
Bi2212^. Two prominent features are observed (Fig. la). 
First, there is a broad maximum around q = (tt, tt). Sec- 
ond, there are incommensurate peaks along the bond and 
diagonal directions which correspond to scattering be- 
tween the antinodal sections of the Fermi surface near the 
(tt, 0) points. As can be seen from the 2D plot in Fig. 2a, 
this scattering forms a square box around the zone cen- 
ter, with maximum intensity at the box corners. These 
peaks in the zero F limit are logarithmically divergences. 
Although potentially important for the charge response, 
they are probably not so for the spin response (since J{q) 
has opposite sign in this region of the zone). 

The other noticeable feature of Fig. la is the broad 
maximum around q — (tt, tt). As a consequence, in the 
gapped case, the dominant behavior below the resonance 
energy is commensurate^^'^^. But this behavior is actu- 
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FIG. 1: Re xo in the normal state at zero energy along the 
symmetry axes of the zone for (a) tbl (b) tb2 (c) tb3 and (d) 
tb4. The zone notation is in tt units. 




FIG. 2: 2D zone plots of Fig. 1. The lower left corner is at 
(0,0) and the upper right corner is at {n,n). The dynamic 
range is the same as in Fig. 1. 



ally the exception, rather than the rule. At this point, a 
discussion on how the fermionic dispersion is chosen is in 
order—. A fit is performed using a six parameter tight 
binding model which includes real space lattice vectors 
of type (0,0), (1,0), (1,1), (2,0), (2,1), and (2,2). In the 
original analysis, the fitting variables were (1) the posi- 
tion of the Fermi surface along (0, 0) — (tt, tt) (the node), 
(2) the Fermi velocity at the node, (3) the Fermi surface 
along (tt, 0) — (tt, tt) (the antinode), (4) the energy at 



FIG. 3: Fermi surfaces corresponding to Fig. 1. The lower 
left corner is at (0,0) and the upper right corner is at (7r,7r). 



(tt, 0), (5) the curvature of the dispersion at (tt, 0) (along 
the kx direction), and (6) the energy at (7r,7r) (this last 
condition is simply for stability of the fit). The behavior 
around q = (tt, tt) in Fig. la is largely a consequence of 
the fact that the saddle point energy at (tt, 0), -34 meV, 
is close to the Fermi energy. For dispersions where this 
is displaced deeper in energy, a minimum is found in Re 
Xo a-t <? = (""j rather than a maximum. 

To illustrate this, a tight binding analysis has been 
performed on the latest high resolution ARPES measure- 
ments on optimal doped Bi2212^^. These data reveal an 
isotropic Fermi velocity with a value of 1.5 eVA. Fixing 
this condition forces the energy at (7r,0) to lie at -119 
meV, as compared to -34 meV in the earlier fit^''. This 
difference (dispersion tb2) has a profound effect on the 
real part of the bubble, as can be seen in Fig. lb. Now, 
incommensurate peaks are observed around q = (tt, tt) as 
well, and as can be seen in Fig. 2b, they form a diamond 
shaped pattern. This structure is actually discernible in 
Fig. la as kinks, but it is obscured by the dominant max- 
imum around q — (tt, tt). The behavior of Fig. lb is also 
found for the t,t' dispersion used in Ref . [25l. as the (7r,0) 
point in that case is at -129 meV. 

The incommensurate effects can also be strengthened 
by flattening the Fermi surface around the nod&i^!^. A 
dispersion which incorporates this effect, as well as the 
presence of a flat quasiparticle band near (tt, 0), has the 
Fermi surface illustrated in Fig. 3c (this is dispersion tb3 
of Ref. m, and is denoted as tb3 here as well). In this 
case, the effect of the enhanced incommensurability and 
the tendency to have a maximum a,t q = (tt, tt) due to the 
small energy of the (tt, 0) point (-10 meV) approximately 
balance, leading to a net flat behavior around q = (tt, tt). 
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as can be seen in Fig. Ic. 

One issue is that there is not much inelastic neutron 
scattering data on Bi2212 because of the small sample 
sizes, though this is starting to change. A large amount 
of data is present for YBCO, but the electronic struc- 
ture of this material is complicated because of its strong 
orthorhombicity and the presence of metallic chain lay- 
ers. Instead, I consider LSCO. This material has the 
advantage of having extensive inelastic neutron scat- 
tering data and a simpler electronic structure (no bi- 
layer splitting). Recently, Tranquada and co-workers^ 
have pointed out that the inelastic magnetic response 
of non-superconducting Lai.875Bao.i25Cu04 (LBCO) is 
very similar to that of superconducting YBCO^^, which 
makes this material of particular interest. Extensive 
ARPES data exist on LSCO. Recently, an in depth study 
was done on underdoped LSCO^S and revealed a quite 
striking Fermi surface (Fig. 3d). The Fermi surface 
is characterized by two straight sections, one centered 
about the node, the other about the antinode. This Fermi 
surface, and the observed dispersion (denoted as tb4), is 
fit using the condition (1) node at (0.44,0.44)7r, (2) ve- 
locity at the node (1.8 eVA), (3) curvature of the Fermi 
surface at the node (zero), (4) antinode at (l,0.18)7r, (5) 
(tt, 0) energy at -70 meV, and (6) (tt, tt) energy at 1 eV 
(for fit stability). The result for the zero frequency real 
response is shown in Figs. Id and 2d. Again, note the 
strong "box-like" structure around the zone center, and 
the pronounced incommensurate behavior around (tt, tt) 
due to the flat Fermi surface near the node. In this con- 
text, it should be noted that the box-like structure near 
q = (0, 0) is due to nesting of the antinodal regions, but 
the diamond-like structure near q = (tt, tt) is due to nest- 
ing of the nodal regions. This is in contrast with the 
alternate stripes model, where the incommensurability 
around these two q vectors have the same origin. 

A summary of the properties of these four tight binding 
dispersions are listed in Table 1. A number of other tight 
binding dispersions have been analyzed as well, including 
some that account for bilayer splitting-^. The results 
fall within the range of behavior discussed for the four 
dispersions here. In addition, bilayer splitting leads to 
some extra details in the even (optic) spin response due to 
differing momentum locations of features from bonding- 
bonding and antibonding-antibonding responses^. For 
simplicity, these details are not discussed here as they do 
not occur in the dominant odd (acoustic) spin channel 
which involves only the bonding-antibonding response. 



IV. RESULTS - SUPERCONDUCTING STATE 

The presence of a d-wave energy gap leads to a po- 
larization bubble which now becomes gapped for all q 
vectors in the first zone but two (these vectors being 
the ones that connect the nodes of the d-wave order 
parameter). Because of the strong anisotropy of the d- 
wave gap, the lower edge of the particle-hole continuum 



TABLE 1: Tight binding dispersions based on angle resolved 
photoemission data. The first four columns list the coeffi- 
cient, Ci, of each term (eV), that is e{k) = ^Cir]i{k). The 
last column lists the basis functions (the lattice constant a is 
set to unity). Below this are several properties of these dis- 
persions - the node, (fejv,fcjv), the antinode, (l,fcAjv) (in vr/a 
units), the nodal velocity, vn, the antinodal velocity, van (in 
eVA), the energy at (vr, 0), em (eV), the bandwidth, W (eV), 
and the value of U needed for resonance at 40 meV in the 
superconducting state. 
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also has a strong q dependence. This continuum edge is 
plotted along the zone symmetry axes in Fig. 4 for the 
four dispersions studied here. The threshold at (tt, tt) 
corresponds to twice the energy gap at the hot spots. 
As one moves away from this wavevector, the threshold 
splits^'^'^^ into several ones: two for q along {tt,x), and 
three for q along {x,x). Along {x,x), the zero threshold 
at qN — '2kN obviously corresponds to two nodal points, 
the analogous minimum along (tt, x) involving only one 
nodal point (the other being at kjsf + (tt, a;)). 

As briefly discussed above, the efi^ect of an energy gap 
has a profound impact on the polarization bubble. The 
presence of an energy gap, along with the fact that the 
coherence factors at threshold do not generally vanish in 
the d-wave case, leads to a step jump in the imaginary 
part of the bubble at threshold. By Kramers-Kronig, 
such a step jump causes a logarithmic divergence in the 
real part of the bubble at threshold. As a consequence, 
any finite value of U will lead to the presence of a col- 
lective mode below threshold. The exception is in those 
cases where and Ak+q have the same sign, which in 
the d-wave case can occur for q vectors significantly dis- 
placed from (tt, tt). In addition, for a pole to occur, then 
U must have a positive sign (this is relevant to near- 
neighbor exchange models, where the positive sign only 
occurs for qx + qy > tt). 

The value of U needed to have a pole at q = (ti", tt) 
at a particular energy varies significantly among the dis- 
persions looked at here (the values for 40 meV are also 
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FIG. 4: Particle-hole continuum edge along the symmetry 
axes of the zone. The maximum d-wave gap is 30 meV. The 
dashed line in (d) corresponds to the process involving fc]v 
and fcjv + q. The zone notation is in tt units. 



listed in Table 1). This is a combination of the different 
Fermi surfaces shown in Fig. 3, the different energies of 
the van Hove singularity at (tt, 0) for these dispersions, 
and the different overall bandwidths. This variation of U 
involves a rather subtle issue concerning whether one uses 
a bare dispersion based on band theory, or a renormal- 
ized (quasiparticle) dispersion, for in Eq. 2 (and also 
the issue of bilayer splitting). Again, the purpose here 
is to simply look at a range of dispersions and contrast 
their behaviors. 

Figs. 5-8 show the variation of the imaginary part of 
the RPA susceptibility (Eq. 1) with U tuned to yield a 
pole at q = (tt, tt) at 40 meV (50 meV for dispersions tb2 
and tb4 - the reason for this difference will be discussed 
below). U was treated as a constant, though calculations 
were also performed for U — J{q). A maximum energy 
gap of 30 meV was assumed. 

Fig. 5 shows results for tbl. At low energies, one sees 
weak intensity near corresponding to the nodes. As 
the energy is raised to exceed the minimum threshold 
along (tTjx), this incommensurate pattern rotates 45 de- 
grees. The preference for this rotated pattern was dis- 
cussed by Schulz^^!^ and is a consequence of the fact 
that translation by (x, x) only brings the Fermi surface 
into coincidence in one quadrant of the zone, but transla- 
tion by (tt, x) brings it into coincidence in two quadrants 
(Fig. 9). In addition, though, one sees a global maxi- 
mum at q — (tt, tt) which is a consequence of the global 
maximum seen in the normal state response at zero en- 
ergy discussed earlier. Above the resonance energy (40 
meV), the response becomes incommensurate, with the 
displacement of the maxima away from (tt, tt) increasing 



FIG. 5: 2D zone plots of Imx for tbl at (a) 10 meV (b) 30 
meV (c) 40 meV and (d) 50 meV. The maximum d-wave gap 
is 30 meV, and U (0.274 eV) is chosen to give a resonance at 
40 meV. In Figs. 5-8, the lower left corner is at (7r/2,7r/2), 
the center at (tt, tt). Scale is such that the lowest is and the 
highest is (a) 0.7 (b) 4.4 (c) 112 and (d) 21.3 



in magnitude (but the response decreasing in strength) 
with increasing energy. For a range of energies above res- 
onance (~ 12 meV), the incommensurate response corre- 
sponds to damped poles, and mostly has maxima along 
along (tt, x), but for some energies the maxima are along 
(x,x). This is due to the anisotropy of the splitting of 
the particle- hole threshold as one moves away from (tt, tt). 
The overall response is best appreciated by plotting the 
intensity as a function of energy along the (x,x) direction 
in momentum space as shown in Fig. 10a. In the pole 
region (40 - 52 meV), one sees a magnon-like (quadratic) 
dispersion. Above this region, the response rapidly loses 
strength and the dispersion becomes more steep. 

I now contrast this behavior with that from the next 
dispersion, tb2. This dispersion, based on ARPES fits to 
the bonding Fermi surface in Bi2212, has (1) an isotropic 
velocity around the Fermi surface and (2) an energy for 
the (tTjO) point which is much deeper than for tbl (-119 
meV as compared to -34 meV). As a consequence, this 
dispersion yields a bubble whose real part has a rela- 
tively weak momentum dependence, and is very typical 
of many of the dispersions used in the theoretical liter- 
ature. In fact, the momentum dependence is so weak, 
that with a constant U, the condition to yield a reso- 
nance at 40 meV at (tt, tt) implies the presence of long 
range order (that is, there is a g value at zero energy 
where the real part of the bubble exceeds l/U). This 
forces us to move the resonance condition close to the 
edge of the continuum (where the real part of the bubble 
has a peak) in order to avoid this problem. And obvi- 
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FIG. 6: 2D zone plots of Imx for tb2 at (a) 10 meV (b) 30 
meV (c) 50 meV and (d) 100 meV. The maximum d-wave gap 
is 30 meV, and U (0.508 eV) is chosen to give a resonance at 
50 meV. Scale is such that the lowest is and the highest is 
(a) 6.3 (b) 20.6 (c) 22.6 and (d) 7.0 




FIG. 7: 2D zone plots of Imx for tb3 at (a) 10 meV (b) 30 
meV (c) 40 meV and (d) 50 meV. The maximum d-wave gap 
is 30 meV, and U (0.214 eV) is chosen to give a resonance at 
40 meV. Scale is such that the lowest is and the highest is 
(a) 4.0 (b) 101 (c) 112 and (d) 10.6 



ously, because of the weak momentum dependence, the 
value of U needed to obtain a resonance condition is sig- 
nificantly larger than for tbl. The results are plotted in 
Figs. 6 and 10b. Below resonance, despite the lack of a 
commensurate response, the incommensurate response is 
similar to tbl, being dominated by node- node processes 
at low energies, and then rotating to a bond centered re- 
sponse once the minimum of the continuum edge along 
(tt, x) is exceeded. The response becomes commensurate 
at resonance (50 meV), then one has a weak incommen- 
surate response above resonance whose maxima generally 
sit along the diagonals. Although this 45 degree rotation 
of the incommensurate maxima is indeed what is seen 
in YBC02^ and LBC021, the incommensurate response 
above resonance still has an overall diamond shape to it 
unlike experiment. This shape is a consequence of 2kF 
scattering which can be seen as well in the normal state 
(Fig. 2b). The overall behavior is best appreciated in 
Fig. 10b, where in contrast to dispersion tbl, a very clear 
'reverse magnon' dispersion is evident for the resonance 
mode. Above resonance, the incommensurate response 
is almost dispersionless, and seems to follow the steep 
particle-hole response involving scattering along the zone 
diagonal (dashed line in Fig. 4d). This steepness is due 
to the high Fermi velocity along the nodal direction and 
unlike the rest of the continuum response, is independent 
of the energy gap. 

The next dispersion, tb3, had been designed to en- 
hance the incommensurability by flattening the Fermi 
surface around the node. In Fig. 7, one can see the strong 
diamond-shaped incommesurate pattern below resonance 



that is a result of the enhanced Fermi surface nesting in 
the nodal region. The evolution of the pattern is again a 
diagonal response at very low energies due to node-node 
processes, a rotation by 45 degrees to a bond centered 
pattern once the threshold along (tt, x) is exceeded, com- 
mensurability at resonance (40 meV), and then again a 
complicated incommensurate response above resonance 
due to the anisotropic splitting of the continuum edge as 
one moves away from (tt, tt). The overall dispersion of 
the magnetic response is most visible in Fig. 10c. Note 
the pronounced downward dispersion of the resonance 
which represents a true pole in the RPA susceptibility 
in this case (this response is below the continuum edge). 
One then sees a weak intensity gap near the node-node 
vector due to crossing into the continuum (the so-called 
silent band'^-). Beyond this, a new pole appears^^ on the 
other side of the continuum (the edge of which is again 
controlled by the high Fermi velocity along the nodal 
direction). One then reenters the continuum, and the re- 
sponse rapidly loses intensity and the damped pole-like 
response is lost above 55 meV. 

I now turn to the last dispersion, tb4. This is based 
on a fit to underdoped LSCO ARPES data^^. The 
Fermi surface in this case is characterized by strong nest- 
ing, both in the nodal region, and also the antinodal 
one. Putting the resonance condition at 40 meV yields 
the same reversed magnon dispersion, silent band effect 
around (node-node vector) , and second mode behav- 
ior for q < qn that was so prominent for tb3 (Fig. 11a). 
The one contrast is at higher energies, the incommensu- 
rability is much better defined than for the other disper- 
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FIG. 8: 2D zone plots of Imx for tb4 at (a) 10 meV (b) 30 
meV (c) 50 meV and (d) 100 meV. The maximum d-wave gap 
is 30 meV, and U (0.276 eV) is chosen to give a resonance at 
50 meV. Scale is such that the lowest is and the highest is 
(a) 1.4 (b) 10.6 (c) 27.1 and (d) 5.5 



(a) 





FIG. 9: Fermi surface for tb4 (dashed line) and its q translated 
image (solid line) for (a) q = (1, .76)7r and (b) q — (.88, .88)7r. 
Lower left corner is at {—it, —tt) and upper right at (7r,7r). 
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FIG. 10: Plots of Imx versus energy and {x, x)tt for (a) tbl 
(b) tb2 (c) tb3 and (d) tb4. Same conditions as Figs. 5-8. The 
intensity is on a logarithmic grid ranging from 0.1 to 100. 
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FIG. 11: Plots of Imx versus energy and {x, x)Tr for tb4 for (a) 
superconducting state with resonance at 40 meV (U=0.332 
eV) and (b) pseudogap state with resonance at 50 meV 
(U=0.343 eV). The intensity is on a logarithmic grid rang- 
ing from 0.1 to 100. 



sions due to the strong nesting (this high energy disper- 
sion follows the steep dispersion along the nodal direction 
shown in Fig. 4d). 

But for most of the results presented here, I choose 
to show instead the case where the resonance condition 
is at 50 meV. This is interesting for two reasons: (1) 
this energy is where the maximum response was seen in 
LBCO^, and (2) this energy corresponds to the contin- 
uum edge at (tt, tt) (and thus where there is a maximum 
in the real part of the bubble). U in this case is adjusted 
so that it is equal to the inverse of this maximum value. 

The results are shown in Fig. 8 and Fig. lOd. In 
this case (as well as when the resonance is at 40 meV), 
the incommensurate response below resonance is much 
cleaner than for the other dispersions, with very clear 



maxima along the (tt, x) direction once the continuum 
edge energy along this direction is exceeded, with x at 
the lowest energies being given by the condition = 
(7r/2 — x/4,7r/2 — x/A). One again sees a clear rotation 
of the pattern by 45 degrees as one crosses through the 
commensurate resonance energy, but as with tb2, the re- 
sponse above resonance has an overall diamond shape. 
Unlike the other dispersions, a well defined incommen- 
surate pattern with maxima along the diagonal direction 
persists to much higher energies, remaining clearly visible 
up to 140 meV. Above this energy, the response is less 
well defined, but still incommensurate. Note the very 
well defined incommensurability at 100 meV, which is 
similar to that of the real part of the bubble (indicat- 
ing pole- like behavior). The interesting point about this 
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FIG. 12; Plots of Imx versus energy and {x,x)tt for tb4 for 
(a) superconducting state and (b) normal state. U is 0.276 
eV. The intensity is on a logarithmic grid ranging from 0.02 
to 20. 



dispersion, as commented above, is that the high energy 
response typically follows the outer branch of the disper- 
sive response corresponding to scattering processes along 
the nodal direction (dashed line in Fig. 4d). This nodal 
velocity is set by the renormalized Fermi velocity along 
the node, which in RVB type models is proportional to 
the superexchange, J. This velocity is also essentially 
constant with doping^, suggesting a connection between 
the universal behavior observed in both ARPES and INS 
data. 

For this dispersion, calculations have also been per- 
formed in the pseudogap approximation. This is shown 
in Fig. lib, and as can be seen, are virtually identical 
to the superconducting results shown in Fig. lOd. This 
indicates that the d-wave phasing relation is not neces- 
sary to obtain the results shown in this paper (just an 
energy gap - though we again note that for the s-wave su- 
perconducting case, there is no resonance effect). In this 
connection, most of the dramatic findings in INS data are 
in the underdoped regime where a pseudogap persists to 
very high temperatures. 

It is also amusing to present results for this dispersion 
when the energy gap is set to zero. This is shown in 
Fig. 12b for the same U as used in the superconduct- 
ing case in Fig. lOd. The response is always incommen- 
surate, and is bond oriented for lu < Uc and diagonal 
oriented for uj > lUc (Fig- 13). lOc is the energy where 
the constant energy contour for the dispersion switches 
from being hole-like to electron-like. Besides the lack of 
a commensurate resonance response, and a spin gap due 
to the energy gap, there are interesting similarities with 
the superconducting results. 



CONCLUSIONS 



o 



FIG. 13: 2D zone plots of Imx for tb4 in the normal state at 
(a) 10 meV (b) 30 meV (c) 50 meV and (d) 100 meV. U is 
0.276 eV. Scale is such that the lowest is and the highest is 
(a) 4.8 (b) 6.8 (c) 7.2 and (d) 3.8 



the particle-hole continuum edge due to the anisotropy of 
the energy gap. There is a striking similarity of many of 
these findings to experimental INS data in LSCO, YBCO, 
and Bi2212 - in particular the reversed magnon behavior 
of the resonance mode, and the finding for several dis- 
persions of a rotation by 45 degrees of the incommensu- 
rate response when passing through the resonance energy, 
which argue for universality in the magnetic response as 
has been commented on in regards to experimental INS 
data. 

On the other hand, there are some differences that need 
to be kept in mind. For instance, the incommensurate re- 
sponse above resonance is generally more variable in the 
calculations than indicated by experiment. The fact that 
this is apparent as well in calculations where the bare 
Greens functions are replaced by experimental Greens 
functions'^^ indicates that this is a general issue. The au- 
thor speculates that the RPA calculations and alternate 
ones based on coupled spin ladders are different limits of 
a more complete theory that properly includes the full 
quantum mechanical nature of both the spin and charge 
degrees of freedom. The development of such a theory 
should help shed more light on the relation of the mag- 
netic fluctuation spectrum and the existence of d-wave 
superconductivity. 



In conclusion, RPA calculations show a remarkably 
rich behavior, in particular in the evolution of the mo- 
mentum response as a function of energy. This behavior 
becomes even richer in the superconducting and pseu- 
dogap phases due to the strong momentum anisotropy of 
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